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ABSTRACT 

We define a large sample of galaxies for use in a study of the fundamental plane in the intermediate redshift 
cluster CL1358+62 at z = 0.33. We have analyzed high resolution spectra for 55 members of the cluster The data 
were acquired with the Low Resolution Imaging Spectrograph on the Keck I 10m telescope. A new algorithm 
for measuring velocity dispersions is presented and used to measure the internal kinematics of the galaxies. This 
algorithm has been tested against the Fourier Fitting method so the data presented here can be compared with 
those measured previously in nearby galaxies. 

We have measured central velocity dispersions suitable for use in a fundamental plane analysis. The data have 
high S/N and the resulting random errors on the dispersions are very low, typically < 5%. Uncertainties due to 
mismatch of the stellar templates has been minimized through several tests and the total systematic error is of 
order -5%. Good seeing enabled us to measure velocity dispersion profiles and rotation curves for most of the 
sample and although a large fraction of the galaxies display a high level of rotation, the gradients of the total 
second moment of the kinematics are all very regular and similar to those in nearby galaxies. We conclude that 
the data therefore can be reliably corrected for aperture size in a manner consistent with nearby galaxy samples. 

Subject headings: galaxies: evolution, galaxies: kinematics, galaxies: elliptical and lenticular, cD, galaxies: 
structure of, galaxies: clusters: individual (CLl 358+62) 



L INTRODUCTION 

Scaling relations have been extensively used to measure the 
evolution of galaxies to redshifts of z 1 . However, without 
mass scales, provided by measurements of rotation curves or 
line-widths, it is impossible to extract the evolution of galac- 
tic M/L ratios from that of the underlying luminosity function. 
The explicit use of galaxy masses, in the fun damental plane 
( [Faber et al. 1987[ iDjorgovski & Davis 1987[ ), or the Tully- 
Fisher relation ( Fully & Fisher 1977 ), makes these scaling rela- 
tions extremely powerful tools for measuring galaxy evolution. 

The fundamental plane (FP) is useful for quantifying the lu- 
minosity evolution of early-type galaxies in distant clusters. It 
is an empirical relation between effective radius, velocity dis- 
persion, and surface brightness, for early-type galaxies. In a 
study of 10 nearby clusters J0rgensen et al. (1996) found 



logr, cx 1.241ogcr-0.821og(/) 



(1) 



with an rms scatter of ±0.084 in logr^ (a 21% distance uncer- 
tainty). The virial theorem and homology imply 



M 

log cx 2 log - log re - log (/) 



(2) 



The implication is that the mass-to-lig ht ratio is a funct ion of 
structural and kinematical observables (Faber et al. 1987), such 
that 

M/Lv cx rO-^VO-^". (3) 



The scatter in Coma of 23% in M/Ly ( F0rgensen ef a/. 1993 ), 
for a given a and r^, suggests a small scatter in age for a given 



(J and re (Prugniel & Simien 1996). It has been postulated that 
systematic departur es from homology exist along th e funda- 
mental plane {e.g., Z!iotti, Lanzoni, & Renzini 1990, but no 



satisfactory model has been created to explain both the slope, 
and the low scatter of the scaling relation. Evolution of the fun- 
damental plane zero-point, slope, and scatter can provide criti- 
cal insight into the evolution of stellar populations of early-type 
galaxies and may shed light on the origin of the fundamental 
plane itself. 

In a pioneering effort, Franx (1995) used fundamental plane 
observations in A665 (z = 0.18) to detect the evolution of stel- 
lar populations in early-type galaxies, van Dokkum & Franx 
(1996) continued this technique with early-type galaxies in 
CL0024+16 (z = 0.39) and measured evolution consistent with 
the passive evolution of stellar populations which formed at 
very high redshift. More recently. Kelson et al. (1997) ex- 
tended fundamental plane measurements to small samples in 
CL1358+62 at z = 0.33, and MS2053-04 at z = 0.58. Their 
results confirmed the evolution in the fundamental plane zero- 
point with redshift. van Dokkum et al. (1998b) have now ex- 
tended fundamental plane measurements to z = 0.83, and have 
constrained both the formation epoch of early-type galaxies and 
Vl. Other authors have also used the elliptical galaxy scaling 
relati ons to constrain galaxy evolution and cosmology as well 



(e. g. , ISchade et al. 1997 



1997, Render et al. 1998 



I 



IHs et al. 199%, Z iegler & Bender 



:^ahre 1998[ ), also finding that E/SOs 
in clusters are uniformly old. 

Both van Dokkum & Franx (1996) and Kelson et al. (1997) 
suggested that the FP tilt may have evolved between z = 
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0.33-0.58 and the present epoch, such that low-mass early- 
types are younger than their high-mass counterparts. Their data 
also indicated that the scatter was compar able to the nearby FP 
(l0rgensenef fl/. 1993, Kelson a/. 1997). However, their sam- 
ples were simply not large enough to accurately determine the 
slope or scatter in the fundamental plane at intermediate red- 
shifts. Kelson et al. (1997) also found that the velocity disper- 
sions of two Eh-A galaxies were low (a ^ 100 km/s). Are such 
low velocity dispersions typical, and are these galaxies rotation- 
ally supporte d, such as wo uld be expected for disk-dominated 
systems {e.g., Franx 1993)? 

In this paper we describe the data and reduction procedures 
used to determine the internal kinematics of a large sample of 
galaxies in the cluster CL1358H-62. The HST data, from which 
the sample was selected, are described in Kelson et al. (1999a) 
and van Dokkum et al. (1998a). Extensive Keck spectroscopy 
has now enabled us to measure velocity dispersions, dispersion 
profiles, and many absorption line rotation curves, for a homo- 
geneous sample of cluster members. In §^, we report on the 
sample selection for the spectroscopy. The observations and 
data reduction procedures are briefly discussed in §|^. The stel- 
lar templates are prepared in §Q. Our method for deriving ve- 
locity dispersions using a direct fitting method is introduced 
in §11 In §||, we test the results against other velocity disper- 
sion routines. In §0, we discuss the strategy for selecting the 
best template and investigate the magnitude of template mis- 
match. In §^ we discuss the aperture corrections required to 
place our velocity dispersions on the same system as the J0r- 
gensen et al. (1996) database, and use the spatially resolved 
kinematics to compare the CL1358H-62 members with nearby 
early-type galaxies. Lastly, sources of error, both random and 
systematic, are discussed in §^. 

2. SAMPLE SELECTION 

We are currently studying, in detail, the galaxy populations 
of three clusters, CL1358H-62 (z = 0.33), MS2053-04 (z = 0.58), 
and MS 1054-03 (z = 0.83), selected fr om the Einstein Medium 
Sensitivity Survey ( Gioia et al. 1992 ). Determining member- 
ship is critical for studying any cluster population. Fabricant et 
al. (1991) and Fisher et al. (1998) have recently collected sev- 
eral hundred redshifts in a 10' x 11' field around CL1358H-62; 
the Fisher et al. (1998) sample is more than 80% complete at 
/?<21mag. 

This field was targeted for extensive imaging with the HST 
WFPC2. A two-color mosaic of the cluster CL1358H-62 cov- 
ering 64n' was constructed using 12 pointings. Of the Fisher 
et al. (1998) database, 194 spectroscopically confirmed clus- 
ter members fall within the field of view of the HST imaging. 
From this large catalog, we selected a sample for detailed study 
with the Low Resolution Imaging Spectrograph (LRIS; O ke et 
al. 1995) at the W.M. Keck Observatory. The high resolution 
of the HST imaging allows us to derive structural parameters 
with the accuracy needed for the fundamental plane (K elson et 
al. 1999a). 

From the catalog of cluster members within the HST mo- 
saic, we randomly selected galaxies down ioR <2\ mag. As is 
discussed in more detail in Kelson et al. (1999a), the /?-band se- 
lected sample contains mostly early-type systems, ranging from 
ellipticals to early-type spirals. Our sample was constructed 
without regard for morphological information, contrary to most 
low redshift samples. The selection was performed with an ef- 
fort to efficiently construct multi-slit plates for LRIS. We used 
three masks, with different position angles on the sky. Dur- 



ing the designing of the multi-slit masks, we utiUzed the appar- 
ent morphologies of several galaxies for major-axis position- 
ing, without impacting the selection criteria. The three mask 
designs overlap in the core of the cluster, and so the sample is 
concentrated towards the cluster center. A total of 67 slitlets 
were allocated for the three masks. 

In Figure |l], we show the color-magnitude diagram for the 
spectroscopically confirmed cluster members within the HST 
mosaic. At the redshift of the cluster, L* corresponds to 7? ~ 
20.1 mag (ado pting the MjLy evolution as measured by K el- 
son et al. 1997). There is clearly a tight color-magnitude re- 
lation in the cluster, as shown by van Dokkum et al. (1998a). 
The entire catalog of spectroscopically confirmed cluster mem- 
bers within the HST mosaic is shown by the circles, and those 
selected for fundamental plane analysis are shown by the filled 
circles. As can be seen in the figure, the fundamental plane 
sample fairly represents th e cluster population, minu s the very 
bluest and faintest objects ( van Dokkum ef a/. 1998a| ). 

Fisher et al. (1998) found 108 members brighter than R = 2\ 
mag within the HST mosaic. Our high-resolution spectroscopic 
sample contains 52 of them. Since three cluster galaxies fainter 
than the magnitude limit were added, the total number of the 
sample presented here is 55. Two of these galaxies were ob- 
served using more than one slit-mask. Eight overlap the sample 
of Kelson et al. (1997). 

This high-resolution sample of cluster members contains 
mostly spectroscopically early-type galaxies, a few Eh-A galax- 
ies, and a few galaxies fainter than R = 2\ mag. The Eh-A galax- 
ies are defined by (H(5+H7+H/3) /3 > 4 A and [Oil] 3727 A < 5 



A (Fisher et al. 1998). The E+A fraction in our high-resolution 



spectroscopic sample is consistent with the fraction Fisher et al. 
(1998) found for the cluster, about 5%. By estimating the mass 
scales for these galaxies, we hope to constrain the relationship 
betw een the Butcher-Oemler effect {e.g.. Butcher & Oemler 
1978, |bressler k (junn l^i^j|,|t^aldwell k Rose 1997i Caldwell 
& R ose 1998) and the histories of " normal" cluster members 
{e.g., [Franx 1993|, [Kelson ef a/. 1997[). 



3. OBSERVATIONS AND DATA REDUCTION 

In May 1996 we obtained high SjN two-dimensional spec- 
tra in the field of CL1358+62. The data were acquired at the 
Keck 1 10m telescope on Mauna Kea using LRIS with the three 
multi-slit aperture plates. The first two were exposed for 8100 
s, and the third for 7400 s. We used the 900 mm"' grating, with 
a dispersion of approximately 0.85 A per pixel. The restframe 
coverage was typically -4100 A to -5200 A. The seeing was 
typicaUy 0'^75 to 0'!8 (FWHM). The slitlet widths were l'^05, 
resulting in a typical resolution of cr, « 1.2 A (-60 kms"'). 
Some slitlets were "tilted" with respect to the primary axis of 
the slit-masks, to align them along galaxy major-axes. These 
slitlets have broader projected widths, and leading to slightly 
poorer spectral resolution. The treatment of the resolution is 
discussed explicitly in §^. 

For the purposes of efficient and accurate extraction of 
one- and two-dimensional spectra, a new set of FORTRAN 
procedures was created. These subroutines are collected in 
a stand-alone package cafled EXPECTOR. ExPECTOR's pri- 
mary functions are that of flat-fielding, rectifying, and wave- 
length calibrating two-dimensional multi-slit spectroscopic 
data. ExPECTOR can be used trivially for long-slit reductions 
as well. 
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3.1. Preprocessing 

Some initial processing is required before execution of 
EXPECTOR. The over-scan and bias-frame subtraction was 
done in IRAF. Cosmic-rays were rejected using median filter- 
ing, with a search algorithm for extended sharp objects. Care 
was taken to not degrade galaxy, or night sky features in the re- 
moval of cosmic-rays. Next, the camera, or "/'-distortion was 
determined by tracing the slit boundaries. This distortion was 
removed from the flat-field images and the cosmic -ray cleaned, 
bias-subtracted data frames. 

The next steps involved flat-fielding the data, which requires 
several stages. First, the spectral dome-flats were normal- 
ized, on a slitlet-by-slitlet basis. The normalization occurs in 
both the X and the y directions and produces two output im- 
ages: a two-dimensional pixel-to-pixel flat-field, and a one- 
dimensional column vector which represents the (multi-slit) slit 
function. After the first pass of the normalization procedures, 
re-normalized columns are medianed to produce a column vec- 
tor which is output as the slit function. The pixel-to-pixel flat- 
field also contains a map of the fringing. For the data presented 
here, however, fringing is negligible. This processing of the 
flat-fields requires that the y-distortion be removed beforehand, 
so that the individual two-dimensional spectra are aligned with 
the dispersion axis of the CCD. 

The data frames are divided by the pixel-to-pixel flat-fields. 
Tests using spectra extending to ~ 9500 A showed that this 
step can be effective in removing fringing. Next, the one- 
dimensional slit functions are shifted (in the y direction), such 
that the slit edges match those in the data frames. These shifts 
are required if the slit edges in the CCD images of the flat-fields 
are not registered with the slit edges in the data frames. Finally, 
the columns in the data frames are divided by the shifted slit 
functions. Each mask has its own set of flat-fields, in order to 
best match the spectral coverage of each pixel, and to match the 
individual slit functions. 

3.2. Rectification and Wavelength Calibration 

The two-dimensional multi-slit images were rectified and 
wavelength calibrated by ExPECTOR using the night sky lines 
([O I], Na I, and OH). A single 2048 x 2048 multi-slit image 
can be fully rectified, wavelength calibrated, and rebinned, for 
example, in just a few minutes on a DEC Alpha processor We 
summarize the procedures below. The interested reader can find 
more details in Kelson (1998). 

The technique for rectifying two-dimensional spectra uses 
Fast Fourier Transforms and efficient cross-correlations to mea- 
sure the pixel shifts between regions of sky (or lamp) lines from 
image row to image row. Thus, our method does not rely on 
centering algorithms to fix the positions of individual lamp, or 
night sky, emission lines. The shifts, in pixels, between adja- 
cent CCD rows are measured as a function of wavelength, by 
performing the cross-correlations in subsections of each spec- 
trum, divided along the dispersion direction. We fit a two- 
dimensional polynomial to the shifts as a function of {x,y) posi- 
tion, and use this polynomial to align all the rows to a common 
coordinate (wavelength) system. The typical rms scatter about 
the fit for these rectification transformations is a few millipixels, 
with several tens of millipixels in the worst (tilted) cases. This 
method for rectifying two-dimensional spectra is extremely ac- 
curate and computationally economical due to its maximal use 
of the available data. 

Because the rectification aligns all the CCD rows, within a 



given spectrum, one only needs to determine the wavelength 
solution for a single row (or average of several, in order to im- 
prove the S/N). Expector's automated wavelength calibra- 
tion routine can use LRIS lamp lines, or the night sky emis- 
sion lines. The algorithm quickly finds the emission peaks, 
measures the FWHM of each peak non-parametric ally, and de- 
fines the center of each line to be the average of the two half- 
power points. When calibrating spectra using night-sky emis- 
sion lines, the algorithm quickly finds an approximate central 
wavelength using cross-correlations, and then identifies triplets 
of PI and P2 OH lines using ratios of spacings between sky 
lines. After an initial determination of a low-order wavelength 
solution, more lines are identified and used to refine the solu- 
tion. 

We chose to fit a fifth order Legendre polynomial for the dis- 
persion solutions. The rms scatter about the fit was typically 
±0.01-0.03 A, though somewhat worse for tilted slitlets due to 
the increase in the width of the lines. Comparison lamp spec- 
tra were used to test the validity of the wavelength solution as 
derived from the night sky emission. The sky lines had equiva- 
lently low rms scatter about the lamp solutions, once the zero- 
point in the wavelength calibration had been corrected. We 
opted for a dispersion solution derived directly from the night 
sky, using emission lines down to 5224 A. 

Each slitlet was rebinned logarithmically to a common wave- 
length range, by chaining the rectification transformations and 
dispersion solutions together. This step involves a single inter- 
polation in the x (wavelength) direction. The output images are 
extra wide to accommodate the staggering of the slitlets in the 
dispersion direction, and flux was conserved in the rebinning 
of each image. Note: the data were rebinned only once in the 
y direction of the CCD and once in the x direction. The expo- 
sures were then averaged, weighting inversely by the expected 
noise, including the noise in the sky and galaxies, as well as the 
read noise in the amplifier. Background subtraction was carried 
out in a manner similar to long-slit reductions, in which the 
sky spectrum at the location of a galaxy is interpolated using 
the observed sky on either side. Residuals around strong night 
sky lines, such as 5577 A, and 6300 A, were interpolated over, 
though this step was performed largely for cosmetic reasons. 

Spectra for the galaxies were extracted by summing five 
CCD rows. The mean S/N ranges from 10 to 80 per A, with 
a median of -35 per A. This high S/N allowed us to sum the 
image rows directly, rather than optimally, which would com- 
plicate the determination of the effective aperture of the ex- 
traction. The spectra themselves are shown in Figure H The 
5 row aperture, 1"08 x 1^^05, is equivalent to a circular aper- 
ture of diameter 1''23 (F0rgensen et al. 1995). In an = 0.5, 
Hq = 65 kms"'Mpc"' universe, this corresponds to a circular 
aperture of diameter -5.5 kpc (or equivalently 11" at the dis- 
tance of Coma). Slitlets tilted by 45° have an effective aperture 
19% larger in diameter 

4. preparation of template spectra 

Because these galaxies are at cosmological distances, one 
must process stellar template spectra with care in order to en- 
sure that they have the the same restfram e instrumental resolu- 
tion as the galaxy spectra (Franx 1993). For nearby galaxies, 
this is obviously not a problem, as the same instrumental setup 
can be used to obtain both the stellar and galactic spectra. How- 
ever, for high-z objects, template spectra cannot be observed 
with the same instrumental setup and one must carefully model 
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the resolution of the galaxy spectra. The stellar spectra must be 
treated with extra caution as wel l if their reso lution is similar 
to that in t he galaxy spectra (see |Franx 1993i v an Dokkum & 
Franx 1996). The spectra of template stars must be convolved 
by a kernel to match the point-spread function of the galaxy 
spectra. 

Using LRIS during twilight in May and August 1996, we 
performed long-slit spectroscopy of several late-type stars (G 
and K stars), listed in Table |l| The 1200 mm"' grating was 
used, providing spectral coverage from 4000 A to 5300 A, with 
a dispersion of approximately 0.64 A per pixel. These resulted 
in spectra with higher resolution that that of our galaxy spectra 
(see § 4-. 3). The telescope was moved in a raster pattern to ex- 
pose the star across the full width of the slit at several spatial 
positions. It is important to fill the width of the slit (0''7), such 
that the spectral resolution is limited by the slit width, and not 
by the seeing. These data were also bias-subtracted, flat-fielded, 
rectified, and wavelength-calibrated. The extracted stellar spec- 
tra have very high S/N, several hundred per pixel. 

4.1. The Wavelength Calibration of the Stellar Spectra 

The stellar spectra were wavelength calibrated using a sparse 
sample of comparison lamp lines in the blue (mostly Ar and 
Hg). A dearth of bright lamp lines in the blue led us to test the 
validity of the dispersion solution when applied to the stellar 
long-slit spectra. We performed a cross-correlation of a high 
resolution solar spectrum, with the twilight and template spec- 
tra, over small, 100 A bins, to provide an independent check 
on the wavelength calibration. There was a high-order system- 
atic trend with wavelength, such that there was a net velocity 
difference between the red and blue sides of the stellar spec- 
tra. The wavelength calibration derived from the internal lamps 
was in error, by about 30 kms"' peak-to-peak, from blue to red, 
in the May stellar spectra, and about 40 kms"' peak-to-peak in 
the August stellar spectra — equivalent to a non-linear stretch- 
ing of about 1 pixel. This can be attributed to the dearth of 
available lamp lines in the blue, coupled with centering errors 
in asymmetric arc profiles, especially in the August lamp data. 
Using this cross-correlation information, we rebinned the raw 
template spectra to the proper dispersion solution. 

4.2. Spectral Resolution of the Template Stars 

We determined the instrumental broadening in our stellar 
spectra in two complementary ways. Initially, the resolution 
of these stellar spectra was determined using the widths of the 
comparison lamp lines. For the May stellar spectra, the a 
widths of these lines are shown in Figure || as crosses. How- 
ever, the internal lamps in LRIS do not necessarily illuminate 
the spectrograph in the same manner as the sky. Thus, we 
wished to compare the arc Une widths with the resolution of 
the stellar spectra themselves. We were able to measure the in- 
trinsic broadening in the stellar spectra by comparing them with 
a high-resolution solar spectrum. 

Using the F ourier Fitting Method (F ranx, Illingworth, & 
Heckman 1989), we determined the broadening required to 
match the high-resolution solar spectrum to the spectrum of the 
twilight sky. The typical resolution in the May 1996 stellar 
spectra was -50 km/s. These Fourier Fitting Method results are 
shown as open circles in Figure ^ for direct comparison with 
the (J widths of the lamp lines. The two methods give similar 
broadening profiles as a function of wavelength. This compar- 
ison gives us confidence that the resolution is well-determined. 



We also used the high-resolution solar spectrum to find the 
instrumental broadening in the stellar spectra themselves and 
found consistent results, except around the Balmer lines, where 
the solar spectrum is clearly not a good match to the later-type 
stars. 

The August 1996 stellar template spectra were taken after an 
upgrade of LRIS, in which a new CCD, field-flattener, and de- 
war were installed. The lamp spectra taken at this time have 
strongly asymmetric profile shapes, probably due to the inter- 
nal illumination ov erfilling the grating ( Gillingham 1996, pri- 
vate communcation). Thus, the spectrograph's resolution of the 
sky (and galaxies) may not be well-described by the arc line 
profiles. Therefore, we used the high resolution solar spectrum 
again, to directly measure the instrumental broadening in the 
August stellar spectra. There was no twilight spectrum from 
August, but the spectrum of the G5IV template star is close in 
type to the solar spectrum, thus minimizing any problems with 
template mismatch in the determination of the resolution. At 
4500 A, the high-resolution spectrum implies an instrumental 
broadening in the HD5268 spectrum of 45 kms"', while the 
lamp line widths yield 55 kms"'. 

If the resolution derived from the August 1996 arc spectra 
were to have been adopted instead of that derived from the so- 
lar spectrum comparison, we would have incurred a -3% error 
in velocity dispersion at <t = 100 kms"'. At cr = 300 kms"', this 
error is only -0.3%. Thus, the effect on the slope, a, of the 
fundamental plane, logr,. cx a log cr + /3 log (/) , would have been 
to increase the coefficient of logcr, by Aa sa +0.02. 

For both the May and August stellar spectra, we adopted the 
resolution determined from the high resolution solar spectrum. 

4.3. The Instrumental Resolution of the Galaxies 

The night sky lines enabled us to measure the variation in 
instrumental broadening as a function of wavelength, for each 
of the CLl 358+62 slitlets. However, care must be taken since 
blended sky lines can bias the measurement of the instrumental 
broadening. Therefore, the spectrograph resolution was mea- 
sured using a list of bright night-sky lines, chosen such that 
they are unblended at our resolution. In particular, we used the 
emission lines at 5577 A, 6300 A, 6364 A, and an array of OH 
PI li nes, all of which are listed in Table ^ (see O sterbrock et al. 
1996). The OH PI lines are close blended doublets, which, for 
our instrumental setup, impacts our estimate of the resolution 
at levels smaller than < 0.5%. The effect is negligibly small for 
our velocity dispersion measurements. We also note that the 
profiles of these sky lines were well fit by pure Gaussians. 

The sky line widths vary as a function of wavelength and 
position across the CCD. We fit third order polynomials to the 
sky-line widths as a function of wavelength, for each slitlet. The 
variation in template resolution with wavelength was also mod- 
eled with a third order polynomial. The final, prepared template 
spectra were generated by convolving the stellar spectra using 
a Gaussian with a a width determined by the quadrature differ- 
ence between the two resolution polynomials. This procedure 
was tested by measuring the effective resolution of the prepared 
template spectra using the high-resolution solar spectrum. The 
effective resolution of the prepared templates agreed with the 
sky line widths to within a few percent. 

Figure ^ shows sky line widths as a function of wavelength 
for each of the CL1358+62 galaxy slitlets in the second mask. 
The circles show the measured widths of the sky lines in kms"', 
while the sohd lines show the polynomial fit to these data. The 
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dashed and dotted lines show the polynomial representation of 
the template spectral resolution for the May and August stel- 
lar spectra (see, e.g.. Figure ||). The resolution for each slitlet 
was determined from the combined spectra of three exposures 
per mask. The widths of the sky lines were also measured in 
the three individual exposures separately to test the stability of 
the resolution. These individual measurements agree with the 
data in Figure ^ to a high degree of accuracy (± 1 %) and are 
therefore not shown. 

Each slitlet was given its own suite of prepared template 
spectra, broadened to the have the same restframe resolution as 
that of the galaxy (or galaxies) contained within. In all cases, 
the resolution of the template spectra was higher than that for 
the galaxies. 

5. VELOCITY DISPERSION FITTING METHODS 

We are endeavoring to find the line-of-sight velocity distribu- 
tion (LOSVD) which gives the best match between our galaxy 
spectrum and a stel lar template broadened by the LOSVD (e.g., 
Rix & White 1992|). We have parameterized the LOSVD by a 



Gaussian, as has been done for comparable studies of nearby 
samples of early-type galaxies. While sometimes the line pro- 
files of real galaxies can be sign ificantly non-Gaussian (see, 
e.g., I^ranx & Illingworth 1988 ), deviations from Gaussian 
LOSVD s are more typically of the order of 10% (B ender et 
al. 1994), and are unimportant for studies of the fundamental 
plane. We use two methods to determine the velocity disper- 
sions, the widely used Fourier Fitting method and a direct fit- 
ting method which is more appropriate for high redshift galax- 
ies (see below). We use both algorithms in order to ensure that 
our measurements are consistent with those in the nearby sam- 
ples. 

5.1. Fourier Fitting Method 

The parameterized LOSVD, has two parameters, the radial 
velocity V, and the dispersion, a, in this velocity. We therefore 
build a model, M, of each galaxy spectrum, G, by a convolution 
M = B(V, a) o T, where T is the template stellar spectrums and B 
is the broadening function. We can then compute the residuals 
between the convolved model and the galaxy spectrum: 



X=\G-M\- 



(4) 



In th e Fourier Fitting method (F ranx, Illingworth, & Heckman 
1989), this calculation is performed in the Fourier domain, 
such that the Fourier transform of the convolved template is 
matched to the transform of the galaxy spectrum: 



X=\G-M[ 



(5) 



When measuring velocity dispersions, we are most inter- 
ested in finding the broadening function which produces the 
best match to the absorption features. While single stellar spec- 
tra can provide good matches to galaxy absorption features, 
the continuum of a stellar spectrum does not match the con- 
tinuum of a galaxy spectrum well at all. For our high redshift 
spectra, the detector response also alters the continuum shape 
with respect to the templates. The low wavenumbers which de- 
scribe the general shape of galaxy and stellar spectra have little, 
if anything, to do with the computational problem of deriving 
LOSVDs, but can contribute substantially to the summation. 

Therefore, an important aspect of measuring velocity disper- 
sions is the matching of the template continuum to that of the 



galaxy. In the Fourier Fitting method (F ranx, Illingworth, & 
Heckman 1989), low wavenumbers are filtered during the sum- 
mation of x^- This filtering helps minimize mismatch between 
stellar templates and galaxy spectra by essentially using sines 
and cosines as additive functions to match the continuum of the 
template to that of the galaxy: 



0, 
[1 
1 



= <{ [l-cos(7r^)]/2, /low < / < 2/io„; 



/low < /; 

/low < / < 

2/low < /■ 



such that the x summation can now be written as 
= l(G-M) X 



(6) 



(7) 



The weight vector equals zero, at frequencies at or below /o„, 
and unity, at frequencies at or greater than 2/iow The filter 
is a continuous function of wavenumber so that the method is 
insensitive to the precise cutoff, and can be thought of as hav- 
ing an effective frequency cut of / = 1 j/ow Put another way, 
continuum matching by the Fourier Fitting Method consists of 
specifically removing long-wave sines and cosines from both 
star and galaxy. In order for our velocity dispersions to be com- 
parable to those of J0rgensen et al. (1995), we adopt the same 
restframe filter window of A: [100]"' and use /ow = 9. 

5.2. Direct Fitting Methods 

When measuring velocity dispersions of high redshift galax- 
ies additional, non-uniform sources of noise in the data can 
complicate the analysis. At high redshift, galaxy spectra are 
redshifted to wavelengths where bright night-sky emission lines 
are present. The shot-noise in these lines are a strong source of 
local noise in the resulting galaxy spectra. Furthermore, sharp 
residuals from the subtraction of the flux in these lines adds 
high-frequency noise to the extracted galaxy spectra. There- 
fore, we constructed a direct fitting algorithm, which allows for 
a more complicated treatment of the data. 

5.2.1. Non-uniform Weighting 

In Fourier methods every pixel receives identical weight in 
the fit of the broadened template to the galaxy spectrum. Uni- 
form weighting, however, is not always desired, especially 
when there exist highly localized sources of noise in the data. 
Such noise adds high frequency power to the power spectrum 
of the data, leading to additional uncertainty in the x^ fit- At 
best, one can attempt to globally reduce the effects of such 
nois e in Fourier methods using, e.g., Wiener filtering (S imkin 



1974, [Bender 1990| , |Rix & White 1992[ponzalez 1993^ . 

By performing the x' computation in the real domain, how- 
ever, one can assign weights to individual pixels. Complicated 
weighting procedures can be adopted, such as, for example, 
masking of regions of poor sky-subtraction, weighting by the 
inverse of the expected noise, or masking of specific spectral 
features such as extraordinarily strong Balmer absorption lines 
or absorption features which have been filled in by emission. 
This process allows such features to be excluded from the x^ 
calculation, leaving the results less sensitive to problematic fea- 
tures and noisy regions of the data. In performing a weighted 
least-squares fit to galaxy spectra, we modify Equation ^ to in- 
clude a vector of weights, W. 



X^=\(G-M)xWf 



(8) 
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In order to properly treat the non-uniform sources of noise 
in our CL1358+62 data, we developed a variant of the Rix & 
White (1992) algorithm which determines line-of-sight velocity 
distributions (LOSVDs) by computing in the real domain. 
As discussed earlier, we adopt a Gaussian function for the pa- 
rameterized LOSVD. 

When using either the Fourier Fitting method, or a direct- 
fitting algorithm, one performs a least-squares fit of broadened 
templates to the observed galaxy spectra. The only essential 
difference between the two methods as applied here is the abil- 
ity to use non-uniform weighting of the data in the real ^ sum- 
mation. Thus, the results of a fit in the pixel (real) domain 
with uniform weighting should be identical to those found us- 
ing Fourier fitting. This is shown to be the case explicitly in 



5.2.2. Continuum Matching 

As with the Fourier Fitting Method, matching the galaxy con- 
tinuum is an important component of our direct fitting method. 
For galaxies at higher redshift, there is a multiplicative source 
of mismatch between the galaxy and stellar continuum as well 
as the standard, additive one. 

When deriving velocity dispersions of nearby galaxies, one 
can obtain template spectra using an identical instrumental 
setup. In this case, any mismatch between the observed shapes 
of the galaxy and stellar spectra arises because the overall 
shapes of galaxy spectra are not perfectly represented using a 
single stellar spectral type; galaxy spectra are complicated com- 
posites of many individual stellar spectra. Such differences in 
the continuum can be approximated by adding a function of 
wav elength comprised o f simple functions, such as polynomi- 
als (Rix & White 1992), to make up differences between the 



flux in the template and the flux in the galaxy. 

The second type of mismatch between the shapes of galaxy 
and stellar spectra occurs because the spectra of distant galax- 
ies are obtained with a completely different instrumental setup 
from the stellar template spectra. In this case, the galaxy and 
stellar spectra may have been obtained using different gratings 
and detectors. Thus, the spectra of the stars and the distant 
galaxies may have different overall shapes because of the in- 
strumental response function, which is multiplicative. 

Therefore, we incorporate a multiplicative polynomial, Pm, 
into the fit to remove large-scale shape differences between the 
observed stellar and galactic spectra which are largely due to 
changes in instrumental setup, i.e., differences in instrumental 
throughput as a function of wavelength. The inclusion of Pu 
in the least-squares solution ensures that our results are insensi- 
tive to the normalization or flux calibration of galaxy and stellar 
template spectra (see below). 

The incorporation of a multiplicative polynomial in the fit 
is a key difference between previous procedures {e.g., Rix & 
White 1992) and this particular direct-fitting algorithm. Instead 
of normalizing the stellar template and galaxy spectra before 
performing the fitting of the template to the galaxy, the nor- 
malization is explicitly incorporated into the continuum match- 
ing process. Since this normalization becomes a part of the 
minimization, the effects of the continuum normalization on the 
derived velocity dispersions are, by definition, minimized. The 
zeroth order term in this multiplicative polynomial is equivalent 
to the "line strength" parameter, 7, in more traditional velocity 
dispersion algorithms. 



The final form of the computation of x is now expressed as 

K 

x'=|{G-[PM(Bor) + P^ + ^fl,//,]}xW|2. (9) 

7=0 

Here, we have written the galaxy spectrum as G and the 
broadened, velocity shifted, template spectrum is written as 
BoT . Once again, we explicitly include W, the vector of pixel 
weights. The coefficients of the M-order Legendre polynomial, 
Pm, are solved for simultaneously along with Pj^ and aj. The 
additive continuum functions can be a polynomial (fV) of order 
A^, or a collection of sines and cosines {Hj) up to order K. 

A simple gradient-search algorithm is employed to search 
for the values of a and V which minimize x^, which includes 
weighting by the inverse of the expected noise. One begins 
with an approximate velocity dispersion, a, and radial veloc- 
ity, y and uses these to broaden and shift the stellar spectrum, 
r. This convolution and velocity shift can be performed in the 
Fourier domain. For a given a and V , it is trivial to compute 
Pm, Pn, and aj, taking into account the pixel weighting. In this 
algorithm, the galaxy spectrum remains unfiltered in any way 
to preserve its noise characteristics and is used directly in the 
computation of x^- When the x^ minimum is found, the errors 
in dispersion and velocity are determined from the local topol- 
ogy of x^(ct,V). 

We use M = 5 to remove the multiplicative shape differences 
between the template and galaxy spectra. Using M = 4 orM = 6 
has almost no effect on the resulting velocity dispersions, less 
than a percent on average. In tests with M = 0, we found that 
the resulting velocity dispersions were systematically higher, 
but by :;5%. When the multiplicative polynomial is included 
in the fit, the reduced x^ values were 25% lower than when ad- 
ditive functions were used alone. This dramatic effect is due to 
the fact that our template spectra have a very strong slope from 
blue to red as a result of the grating response function. Match- 
ing the overall shapes of the templates with those of the galaxy 
spectra is clearly an important step in the procedures. For work 
on distant galaxies, the difference between the instrumental re- 
sponse function of the galaxy spectra and the templates should 
not be removed using additive functions alone. By incorpo- 
rating the renormalization into velocity dispersion fitting pro- 
cedure, we can be confident that all instrumental effects have 
been accounted for, and their effects on the results have been 
minimized in a least-squares sense. Given the above test, we 
are confident that differences in the instrumental setup lead to 
systematic errors which are ^5%, and are probably of order 
-1%. 

The parameter K is equivalent to the /io„ parameter of the 
Fourier Fitting Method when one uses sines and cosines as the 
basis functions for the additive continuum components. As dis- 
cussed earlier, /io„ is the half-power point of the filter in the 
Fourier domain, and the effective number of frequencies re- 
moved from the x^ sum is K = Ij/iow For these CL1358H-62 
data, we use K = \1>. For consistency with this filtering, the 
number of polynomials one must use in the continuum fitting 
is defined by the number of zeroes in the sine and cosine fil- 
ter. Since each sine or cosine of frequency / has 2/ zeroes, 
N = 3/iow The results are insensitive to the choice of at the 
level of a few percent, and merely reflect our susceptibility to 
template (continuum) mismatch, which is discussed below. 



5.2.3. The Generalized Direct Fitting Procedure 



5.2.4. Expanded Variants of the Direct Fitting Method 
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Additional terms can be added to Equation ^ for a variety of 
specialized cases. For example, multiple template spectra can 
be incorporated straightforwardly: 



x' = |{g-[Pa 



N, 
k 



(10) 



In such a case, the fitting is complicated by the requirement of 
positive definite solutions for all 7^, the relative contribution of 
each template, T^. Each template may have its own broaden- 
ing function B^, such that each component has its own radial 
velocity and velocity dispersion.. 

Furthermore, for data in which subtraction of the sky back- 
ground is a difficult procedure, a high S/N sky spectrum, S, 
may be incorporated into a fit directly to the raw data, D. Thus, 
the raw data are approximated by a a sum of the sky with the 
broadened template spectrum: 



K 

= \{D-[PLS+PM(BoT) + PN+Y,^jHj]} X W\\ (11) 



In this case, instead of a sky-subtraction being performed on an 
unconstrained column by column basis, as is nominally done, 
each pixel in the spectrum contributes to the fit. The sky spec- 
trum is scaled by a low-order polynomial Pl of order L. In the 
simplest case, Pl can be a constant, appropriate when large- 
scale flat-field variations have been accurately removed. How- 
ever, when large-scale flat-fielding may be an issue, the order L 
should be increased. One can also fit for small shifts in the the 
sky vector if small rectification errors remain in the data. Such 
an application may prove useful when working with fiber-fed 
spectrographs, where determination of an accurate local sky 
spectrum is problematic, or when analyzing two-dimensional 
spectra of the l ow surface brightness halos of cD galaxies (K el- 
son etal. 1999d). 

6. TESTS BETWEEN FITTING METHODS 

Our new direct fitting method has been tested against another 
direct fitting algorithm, a direct analog of the Franx et al. (1989) 
Fourier Fitting Method which was used by van Dokkum et al. 
(1998b) to measure velocity dispersions for galaxies in the clus- 
ter MS 1054-03. They used a program analogous to the Fourier 
Fitting Method, though its treatment of the continuum filtering 
is not given by the smooth cosine bell, but has a rigid cutoff 
instead. This treatment is similar to the Rix & White (1992) 
prescription for treatment of the continuum for the case where 
the continuum functions are sines and cosines. 

Since our template spectra had different shapes than our 
galaxy spectra, due to the overall efficiency of the spectrograph 
decreasing to the blue, we were required to re-normalize the 
spectrum of the stellar template to the overall shapes of the 
galaxy spectra. The new methodology which we described in 
§^ already includes a multiplicative polynomial in the fitting 
process. The Rix & White methodology does not, nor does the 
Franx real-fitting variant of the Fourier Fitting method since 
these programs only use additive continuum functions in the 
fitting procedures. 

We have compared our new direct fitting algorithm to the real 
variant of the Fourier Fitting method. The aperture-corrected 
values (§||) are listed in Table ^ in which the results from the 
method of § 5^ are referred to using ctdf.k and values derived 



using a real-fitting variant of the Fourier Fitting method are re- 
ferred to by (Jdf.f- After renormalizing the templates to the 
shapes of the galaxy spectra, and deriving new values of a with 
Franx's real-fitting variant of the Fourier Fitting method, we 
found a mean (median) offset of -0.6% (-0.4%) compared to 
the results of the new direct fitting algorithm described above, 
with a standard deviation of 2%. This comparison is illustrated 
in Figure |], in which we show a histogram of the difference 
between logcr using the two codes. Another very important 
consistency check is the comparison of the reported formal er- 
rors, da- The two programs report the same formal errors in the 
mean to < 1% of the error, with a standard deviation of 5% of 
the error 

In order to compare results from our new direct fitting 
method with that of the original Fourier Fitting Method, we 
turned off all pixel masking in the direct-fitting procedure. 
Thus, any strong residuals due to sky-subtraction needed to be 
interpolated over to ensure that they did not strongly affect the 
velocity dispersion fitting, corresponding to the usual proce- 
dures for Fourier fitting. 

For templates with spectral types later than solar, agreement 
is quite good, with very small systematic differences between 
velocity dispersions measured with the two algorithms. For 
most of the stars, the agreement is at the level of 1 %, with 
a typical scatter of 3-4%. Such offsets can arise from differ- 
ences in the way in which the stellar continuum is filtered and 
matched to the galaxy spectra. The continuum filtering is an 
important part of minimizing the mismatch between the broad- 
ened stellar templates and the galaxy spectra, and any differ- 
ences in this process will be most important for stars which 
already provide poor matches to the galaxy spectra. In Figure^ 
we show the comparison using the adopted template HD72324, 
for those galaxies which required no masking of emission or 
excess Balmer absorption. 

After these tests, we conclude that the velocity dispersion 
measurements made using the direct fitting method are com- 
parable to those made with both the classical Fourier Fitting 
Method and its real-fitting variant, as used by van Dokkum et 
al. (1998b). 

7. MINIMIZING TEMPLATE MISMATCH 
7.1. Selecting the Best Stellar Template 

The fitting procedure(s) outlined in the previous sections pro- 
vide two diagnostics for determining which template star best 
matches the galaxy data. First, the direct fitting procedure pro- 
duces a fitted template, which allows one to visually inspect the 
residuals and the quality of the fit. In general, visual inspec- 
tion allows one to sort out the best two or three templates. The 
determinations allow us to determine the best-fit template 
quantitatively, rather than visually. 

Many of the galaxies show enhanced Balmer absorption, two 
show emission, and sharp sky-line residuals exist in nearly all of 
the galaxy spectra {e.g. 5577 A). Therefore, for every galaxy, 
we chose to consistently mask the regions around Na and [O 
I], the strongest sky-line residuals in the spectral regions fit- 
ted. Some galaxies also required the masking of strong Balmer 
lines. The Eh-As {e.g. ID# 209, 328, and 343) have strong, 
broad H/? and H7 absorption which are not well-matched by 
any of our current set of template spectra. The emission line Sb 
galaxy ID# 234 and the cD galaxy ID# 375 required the mask- 
ing of the H/3 and H7 emission features, as well as [O III]. 

The stai- with the lowest median x^> overall, is HD72324 
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(G9III). We therefore adopt the velocity dispersions found 
using the template HD72324, regardless of any individual 
galaxy's lowest -y^ . The star which gave the second low- 
est was HD102494 (G9IV). Dispersions derived using this 
star were only offset systematically from values derived using 
HD72324 by 1%. 

Below each galaxy spectrum in Figure ^ we plot the residu- 
als from the fit of the template HD72324. Regions which have 
been masked because of sky line residuals are also masked from 
these plots. Regions of the spectra which were excluded from 
the fitting due to emission or abnormally strong Balmer absorp- 
tion are shown in the figure, and the reader can easily see that 
the galaxies possess a wide range of stellar population and star- 
formation properties. The Balmer and metal line characteristics 
of the samp le will be analyzed in detail in a later paper ( Kelson 
etal 1999c). 

There still remains the possibility of residual template mis- 
match within the sample, artificially inflating the scatter in the 
fundamental plane. To test the magnitude of this effect, we 
compared dispersions derived using HD72324, to the set of dis- 
persions found using the lowest on an individual basis. For 
those galaxies which do not have HD72324 as best-fit template, 
the mean offset is 1%, with a standard deviation of 2.8%. Since 
these galaxies make up less than half the sample, the effect is 
actually much smaller over the entire sample. 

7.2. Template Mismatch 

Our method of determining velocity dispersions explicitly 
assumes that a single stellar spectrum can provide a good match 
to any given galaxy spectra. In general, one should not expect 
any single star to give a perfect fit, since galaxy spectra are com- 
prised of the spectra of a range of stellar spectral types. Some 
of the mismatch between the stellar template and the galaxy 
spectrum is compensated by the continuum matching described 
earlier Remaining mismatch over spatial frequencies smaller 
than the filter window, both in the continuum, and in the de- 
tailed depth and shapes of the absorption lines, can be fitted 
using com posite templates. Wh ile the fit can be improved sig- 
nificantly ( Rix & White 1992 ), this effect is, at most, a few 
percent ( Gonzalez 1993). 

Discrepant absorption line depths and small-scale power in 
the continuum will depend on several physical parameters, the 
overall chemical and spectral makeup of a galaxy itself, the 
metal abundances in the star from which the template spec- 
trum was obtained, the spectral type of the star, etc. Therefore, 
one might expect that the velocity dispersions one finds will 
vary, systematically, with the spectral type an d metallicity of 
the adopted template (see, e.g., Gonzalez 1993). 



7.2.1. SpectralType 

Figure |^ shows the sensitivity of the measured velocity dis- 
persions for our galaxies to the template's spectral type in the 
direct-fitting method. The figure shows how the variation of a 
with template spectral type is systematic. The separate panels 
show the logarithmic difference between each of the templates 
and the average of the five G9-K3 stars, chosen as a reference. 
The most deviant template is that of the solar spectrum (G2V). 
Using the solar spectrum as a template gives dispersions which 
are quite deviant, 15% lower on average, systematically. The 
next worse template is the G5IV star HD5268, which gives 
values 4% lower than those obtained with the best fit template 
HD72324. 



A variation of mismatch with dispersion can be a worry for 
fundamental plane work, since, in principle, the slope with re- 
spect to logo- can be biased. In the extreme case of the solar 
spectrum, the offset in a from the values obtained from late- 
type stars differs between a lOOkms"' galaxy and a 300 kms"' 
galaxy by about Alog(o') « 0.1; the slope of the FP is then bi- 
ased by 5a w 0.2. However, the stars with the lowest values 
yield values for a which systematically differ by at most a per- 
cent or two over such a long baseline of a. Thus, the impact on 
the fundamental plane slope is only a few percent. 

7.2.2. Metal Abundance 

Template mismatch is also a function of the absorption line 
strengths in the stellar spectra (at a given spectral type). Two 
of our stars, th e two G9III stars, have m etallicities differing by 
a factor of 1 (t^uck & Challener 19951 P ilachowski, Sneden, 
& Kraft 1996). To zeroth-order, this difference is compensated 
for by 7, the line-strength parameter in the velocity dispersion 
fits. This compensations is clearly seen in the case of these 
two stars, as the average ratio of 7hd72324 /7hd6833 ~ 0.8. These 
two stars with identical spectral type, and a factor of 10 dif- 
ference in metal abundances, show a 3% systematic offset in 
velocity dispersion. They also show a 3% systematic trend in 
a over a baseline of 0.5 dex in a, but only a 2% trend when 
one excludes those galaxies with a < 100 kms"'. By choos- 
ing the best-fit template, one's systematic template mismatch 
errors should therefore be smaller still. By simply choosing a 
late-type star, one is already reducing the mismatch errors to a 
few percent. 

Procedures have been developed which combine multiple 
template spectra to c reate an optimal match to the galaxies (e.g., 
Rix & White 1992). Such algorithms adjust the mix (i.e., 7) 



of stars to create a non-negative distribution of stellar compo- 
nents. The direct-fitting method can be expanded fit foran op- 
timal combination of template spectra, as shown in §5.2.4 but 
such an algorithm is not necessary for our purposes. Our uncer- 
tainties due to template mismatch are already small, at a level 
comparable to the formal errors from the velocity dispersion 
fits. 

7.3. Wavelength-Dependent Mismatch 

There is an additional important test which serves two pur- 
poses. By comparing velocity dispersions derived from dif- 
ferent sections of the galaxy spectra one can determine if (1) 
one has indeed used a template which best matches the galaxy 
spectra globally; and if (2) the measured velocity dispersions 
depend on the spectral region used in the fit. The latter is an 
important question since our fitting is targeted at a region of the 
spectrum containing the G band; we wish to know if our mea- 
surements can be compared to velocity dispersions measured in 
other regions of galaxy spectra (e.g., Mg b). 

If one is using a template which provides a perfect match 
to a given galaxy spectrum, then one should measure identical 
values for a using any region of the spectrum. The practical ap- 
proach to this is to test if the dispersion measured from the blue 
side of the spectrum equals that measured from the red side. 
Therefore, we split our galaxy spectra in two, and measured a 
from the two halves, for the 46 galaxy spectra with S/N > 10 
per pixel in both halves. 

Figure ^ shows the logarithmic difference in a from the two 
sides for each template. The best-fit template star, HD72324, 
yields a mean difference in a between the two halves of less 
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than a percent. In fact, most of the template stars give reason- 
able results. The rest of the stars, aside from the solar spectrum, 
all show absolute differences on the order or less than a few per- 
cent. The solar spectrum shows a median difference between 
the two halves of -10%. 

As a result of these multiple tests, we are confident that our 
errors due to template mismatch are likely to be small, at the 
level of a few percent, comparable to the formal errors in the 
fit. 

8. THE VELOCITY DISPERSIONS 

Early-type galaxies have radial gradients in velocity disper- 
sion and radial velocity. The implication is that the a one mea- 
sures depends upon the metric aperture used to extract the orig- 
inal spectra. The aperture is defined by an angular extent upon 
the sky, so the metric aperture which defines the velocity dis- 
persion typically increases with distance to a given galaxy. 

The dependence of a upon aperture is not straightforward 
to predict because the measurement is a luminosity-weighted 
mean velocity dispersion within the predefined aperture. Fur- 
thermore, early-type galaxies can be partially supported by ro- 
tation, and velocity dispersion measurements made using large 
metric apertures will reflect the total second moment of the inte- 
grated line-of-sight velocity distribution: /i2 = V a^ + V^. The 
integration of this quantity is weighted by the surface bright- 
ness distribution of the galaxy within the aperture of the spec- 
trograph. Thus, measurements of a depend on each galaxy's 
intrinsic distribution of orbits, as well as its light distribution. 

8.1. Spatially Resolved Kinematics 

Our velocity dispersion measurements will be compared di- 
rectly with similar measurements made with nearby galaxies, 
e.g. J0rgensen et al. (1996). Therefore, one source of concern 
is that the galaxies in the CL1358H-62 sample may have differ- 
ent internal structure than nearby galaxies. We begin exploring 
the magnitude of this effect by comparing the kinematic profiles 
of our galaxies with those of nearby galaxies. 

With the excellent image quality and high S/N of the data, 
we were able to extract rotation curves and velocity dispersion 
profiles. Using the direct-fitting method of we measured 
V and a along each slitlet. When necessary, rows were binned 
to achieve sufficient S/N. The sum used the same masking 
and noise-weighting as in the determination of the In" aper- 
ture dispersions. Two of the galaxies have S/N ratios too poor 
to provide spatially resolved kinematics (# 360 and # 493). 

The profiles can be seen in Figure ^ where absorption line 
rotation curves and velocity dispersion profiles are plotted ad- 
jacently for each galaxy. These profiles have not been corrected 
for seeing, nor for non-major axis positioning of the slitlets. For 
most high-redshift galaxies, the rise in the rotation curve has 
been heavily modified by seeing (see, fo r example, the model- 
ing of|Vogtef g/. 1994 |Vogt et al. 19971). 



8.2. Comparison with Nearby Galaxies 

In order to test whether aperture corrections derived from 
nearby galaxies can be applied to the high redshift sample, 
we now compare the kinematic profiles in Figure ^ to those 
of nearby galaxies. For this comparison, we use 83 early- 
type galaxies in the kinematic database of Simien & Prugniel 
(1997a,b,c) and construct /i2 = Vcr^ + V^ as a function of ra- 
dius using their long-slit rotation curves and velocity disper- 



sion profiles. In Figure we show the logarithmic gradi- 
ents of /i2, normalized by the mean value within r^/lO (we plot 
only those data with uncertainties smaller than 10%). A least- 
squares fit to the data in Fig. |l^ yields a logarithmic gradient 
of c/ ln(/i2)/i^ log r = -0.023 ±0.006 per dex, shown by the solid 
line. The effective ±1 and ±2(t scatter about the fitted gradient 
is shown by the dashed and dotted lines, respectively. 

In Figure |ri](a,c), we plot the profiles of Vcr^ + V^, con- 
structed from the data in Figure |, for all of the CL1358H-62 
galaxies. As in Figure |lO[ we plot only those data with un- 
certainties smaller than 10%. These gradients have not been 
corrected for seeing. In (b,d) we show the gradients, normal- 
izing the position along the slit by the galaxy half-light radii, 
taken from Kelson et al. (1999a). In these figures we plot the 
first-order fit from Fig. |l^ as the solid line. We show the ±lcr 
(dashed) and ±2cr (dotted) contours of the nearby galaxies as 
well. A least-squares fit to all of the CL1358H-62 galaxies yields 
a logarithmic gradient of (iln(/i2)/iilogr = -0.017 ± 0.026 per 
dex, consistent with the sample of nearby galaxies. 

From the figures, we draw several important conclusions. 
First, the p2 profiles of E/SOs and early-type spirals have small 
gradients, both locally, and at z = 0.33. Second, within r < rif2, 
the gradients are independent of morphology. Least-squares fits 
to the gradients in the spirals are not significantly different from 
fits to the gradients in the E/SOs. The scatter in the p2 profiles at 
z = 0.33 is remarkably consistent with that in the nearby galax- 
ies. We therefore conclude that aperture corrections derived 
from nearby galaxies can safely be used on the full sample of 
CL1358H-62 galaxies. 

8.3. Aperture Corrections 

J0rgensen et al. (1995) used existing spectroscopic and pho- 
tometric data of nearby galaxies in the literature to construct 
models which were "observed" through a series of concen- 
tric apertures of increasing size. Major-axis long-slit slit spec- 
troscopy, by itself, is not suitable for deriving such corrections. 
They found that the following power law can be used to correct 
the observed dispersions, a^, to a physical aperture: 



log (Tc = log(Tz + 0.04(log^^7-log^ic), 



(12) 



where CTc is the corrected velocity dispersions within an aper- 
ture dc, and d^ is the aperture of the observation. For this paper, 
we use dc as the nominal aperture of 3''4 for galaxies at the dis- 
tance of Coma, and d, is the effective aperture for CL1358H-62 
of 1 "23 (see §1). 

We used Eq. |l2| to correct the observed velocity dispersions 
to an aperture consistent with the comparison sample of J0r- 
gensen et al. (1996). Using q^) = 0.5, we obtain values from 
the standard aperture at Coma of dc = 1.68 kpc, and an aper- 
ture size of d, = 5.50 kpc at z = 0.33. Therefore, our measured 
velocity dispersions underestimate the central values by 4.7%. 
Most of the slitlets were not tilted and thus do not suffer addi- 
tional broadening. The most extreme case of a 45° tilt (see §H) 
changes the aperture correction to 5.4%. Thus, we adopt a sin- 
gle aperture correction, regardless of slit tilt. The final velocity 
dispersions have therefore been multiplied by a factor of 1 .047. 

The velocity dispersions, corrected to an aperture equivalent 
to the nominal aperture at the distance of Coma, are listed in 
Table pi A histogram of central velocity dispersions is given in 
Figure]l2|. 
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9. SOURCES OF UNCERTAINTY 

When measuring velocity dispersions potential sources of er- 
ror can exist at several stages in the procedures. Uncertainties 
are introduced by the templates themselves, as well as by the 
galaxy spectra. Poor S/N clearly introduces large random er- 
rors, but can also introduce systematic errors ( J0rgensen et al. 
1995). Furthermore, the imposition of a single stellar compo- 
nent on a galaxy spectrum, of any S/N, can produce mismatch 
errors. We attempt to itemize these sources of error in a quan- 
titative manner, and estimate to what extent our measurements 
are affected by them. 

9.1. The Instrumental Resolution 

At the outset, a potentially large source of error can be im- 
properly measured resolution. In §^ we found that the uncer- 
tainty in the absolute resolution of the template spectra is un- 
certain by no more than ~5 kms"'. This translates to an error at 
a = 100 kms~' of -2.5%. Since the resolution is more important 
for low velocity dispersions, there is almost no effect for galax- 
ies with a = 300 kms"' . Over the baseline of 100 kms"' to 300 
kms"', the mean systematic effect on the velocity dispersions 
is of order 1-2%. Given that the systematic effect is a function 
of velocity dispersion, there is a net uncertainty in the derived 
fundamental plane slopes of about A(9 log r(./9 log cr) k, ±0.02. 

9.2. Aperture Corrections 



In §8.2, we showed that the total kinematic profiles of the 
C11358H-62 galaxies is similar to those of nearby galaxies, and 
concluded that the J0rgensen et al. (1995) form of the aperture 
corrections remained valid for our sample. We therefore con- 
clude that any systematic error in adopting their aperture cor- 
rections for the entire sample is going to be negligible compared 
to other sources of error, such as template mismatch. More so- 
phisticated modeling is required to better estimate this uncer- 
tainty. 

9.3. Signal-to-Noise Considerations 

According to J0rgensen et al. (1995), if the S/N ratio of 
galaxy spectra are too low, the derived velocity dispersions will 
be measured systematically high (see their Figure 3). Five of 
our spectra have S/N i,2Q per A (two are of the same galaxy, 
#493). Only one has S/N < 15 per A. It is #360 with a mea- 
sured (7 = 173 ±20 km/s (a 12% uncertainty). According to 
J0rgensen et al. (1995), this measurement is in error by only 
H-2%. For this galaxy, the low S/N appears to be due to poor 
subtraction of the background in a crowded region of the sky. 
The slitlet containing #360 is 12" long, with #360 at one end, 
partially covering the outskirts of the cD #375 at the other end, 
and containing #368 in between the two. The other four spec- 
tra have measured velocity dispersions between 100 km/s and 
140 km/s and, according to the simulations of J0rgensen et al. 
(1995), may be systematically too large by 1-2% as well. Given 
that these systematic uncertainties are much smaller than the 
formal uncertainties, and that only four galaxies, at most, are 
affected, we do not include any corrections for this effect. 

9.4. Template Mismatch vs. Signal-to-Noise 

The formal errors in the velocity dispersions have two pri- 
mary components: random errors due to photon statistics and 
systematic errors due to, for example, template mismatch. In 
this section, we determine to what extent our formal errors are 
dominated by the random or systematic eiTors. 



We compare the systematic errors due to template mismatch 
to the random errors by comparing three estimators of S/N: (1) 
the S/N ratio per A expected solely from considerations of the 
read noise and noise in the sky and galaxy; (2) <j/Sa, the ratio of 
the velocity dispersion with its formal error; and (3) Q, the ratio 
of the mean flux level in the galaxy and the standard deviation 
of the residuals in the velocity dispersion fit. In the case where 
one uses a perfect template spectrum and no noise has been 
added to the data during the reduction process, Q should be 
equivalent to the expected S/N. All of these are related (since 
the expected noise was used in the calculation of weights in the 
fitting), and we can use them to disentangle the contributors to 
our uncertainties. 

In Figure [T^a) we plot 2 as a function of the expected S/N 
ratio (per pixel in this case). The data do not follow the line of 
unity slope. For those spectra with high S/N, the residuals from 
the fit are larger than would be expected from photon statistics. 
This point is repeated in Figure 13^b) in which we specifically 



show the ratio of Q/ {S/N) as a function of expected S'/A^. When 
the S/N approaches -30, the template mismatch errors become 
the dominant source of uncertainty. Therefore, we do not see a 
constant level in (b), but a trend with S/N in which the quality 
of the fit begins to suffer compared to what one would expect 
from noise considerations alone. This fact can also be seen by 
inspecting the residuals from the fits in Figure ^ in which there 
appear to be consistent features in the galaxy spectra which are 
not well matched by the the fitted template. 

In figure |l3|(c) we plot (t /6a, the formal signal-to-noise in 
the dispersion, as a function of Q, the quality of the fit for each 
spectrum. For most galaxies, there is clearly a one-to-one cor- 
respondence. In (d), the ratio is plotted as a function of Q so 
that one can more clearly see the discrepant data points. Thus 
our error estimates are consistent with the rms scatter about the 
fit, though some galaxies clear deviate from this expectation. 

There are a handful of spectra, however, with velocity disper- 
sion uncertainties which are larger than expected from the rms 
residuals of the fit. We show in in Figure |l3|(e) that the errors 
are dominated by atypically large template mismatch. These are 
the galaxies with small values of the the line-strength parame- 
ter 7. These points with very low 7 values are galaxies #209 
(both spectra), #234, #328, and #343, the Eh-A and emission- 
line galaxies. For these cases, the formal errors in the velocity 
dispersion are higher than expected for the quality of the fit of a 
single template star, even though the strong Balmer absorption 
lines and emission lines were given zero weight in the fit. 

We conclude that for the galaxies with S/N > 30, our veloc- 
ity dispersion uncertainties are dominated by systematic errors, 
such as template mismatch errors (and as noted, these errors are 
^5%). 

9.5. Internal Consistency 

For the galaxies in common with Kelson et al. (1997), the 
agreement in cr is 1% in the mean, with a scatter of 4-5%, con- 
sistent with the formal uncertainties. 

We can estimate the internal consistency of the dataset us- 
ing two galaxies which were measured in two different masks, 
galaxy ID# 493 was observed at two different position angles, 
45° apart. The two velocity dispersions, 117 ±8 kms"' and 
128 ± 9 kms"', agree very well. The other galaxy, ID# 209, is 
an Eh-A whose spectrum required masking of the strong Balmer 
absorption. The two observations for this galaxy do not agree 
because the one with the smaller a was made with a slightly 
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mis-aligned slitlet. In this discrepant observation, the galaxy 
was a secondary target in a slitlet which was centered on a dif- 
ferent galaxy and subsequently tilted to include # 209. The slit- 
let was not correctly tilted to encompass the center of # 209, and 
thus the galaxy was not centered in the slitlet. Thus, the seeing 
convolved galaxy did not evenly fill the slit, leaving a strong 
gradient of its light profile across the slit. Thus, its light pro- 
file under-filled the slit width, leading to an improper estimate 
of the resolution of its spectrum. Furthermore, the resulting 
dispersion is not a central value, but one offset from the center. 
Thus, the larger value of a= 108 ±5 kms"' is the more credible 
measurement. No other galaxies in this sample were secondary 
targets within slitlets. 

9.6. The Effects of Under-filling Slitlets 

At these high redshifts, it might be expected that small galaxy 
sizes, combined with the good seeing of our observations, 
might lead us to under-fill the 1"05 wide slitlets. If a galaxy 
under-fills a slitlet, then the broadening, as measured by the 
sky-line widths, would be overestimated, because the galaxy 
light is coming from an effective aperture narrower than the 
slitlet. Fortunately, the galaxies in this sample are generally 
not small enough to cause problems. As shown in Kelson et 
al. (1999a) the smallest galaxy, #493, has a half-light radius 
of ~0"2. Using the spatial profile along the slit of this galaxy, 
within a box 1"05 wide, we have convolved a spectrum of the 
sky, and compare the widths of the convolved lines with the 
widths of the sky lines convolved by a top-hat 1''05 wide. One 
does find that the sky lines are narrower in the case where the 
galaxy profile was used in the convolution, but the effect is ap- 
proximately 2% of the measured widths. We therefore conclude 
that under-filling of the slitlets is a negligible source of error 

Given the number of sources of systematic error, we estimate 
that the total systematic errors in the order of 5%, generally 
larger than the random errors expected from photon statistics 
alone 

10. SUMMARY AND CONCLUSIONS 

We have accurately measured the internal kinematics of a 
large, homogeneous sample of galaxies, within the cluster 
CL1358H-62. These data will be combined with the structural 
parameters in Kelson et al. (1999a) for the analysis of M/L ra- 
tios and the fundamental plane of early-type galaxies ( Kelson 
et al. 1999b. These direct measurement of velocity dispersions, 
will provide crucial mass estimates for these galaxies, enabling 
us to provide constraints on early-type galaxy evolution and 
their epoch(s) of formation. The formal errors are quite small, 
typically a few percent. For galaxies with S/N > 30 (roughly 
half the sample), the uncertainties are dominated by the system- 
atic errors, which appear to be at a level of 5%. 

These velocity dispersion measurements can be compared di- 
rectly to those made in nearby galaxies. We have measured a 



using a new variant of direct (real) fitting. Results from this 
method have been compared with those of the Fourier fitting 
method ( |^ranx, lllingworth, & Heckman 19*8^ , and we find 
excellent agreement between the two programs. Differences 
between the two algorithms lead to systematic uncertainties in 
(7 on the order of ~±1%, with a scatter consistent with the for- 
mal errors. 

We have measured absorption-line kinematic profiles for 
most of the galaxies in this sample. The CL1358H-62 galaxies 
show a wide range of rotational support. However, the profiles 
of the total second moment of the LOS VDs, /i2 = V cr^ + V^, do 
not vary greatly from galaxy to galaxy. Thus, we conclude that 
a single aperture correction can be applied to our entire sample. 
Furthermore, because the kinematic profiles of the CL1358H-62 
galaxies are similar to those of nearby galaxies, we conclude 
that we have reliably corrected our large-metric-aperture ve- 
locity dispersions to an effective aperture consistent with the 
nearby galaxy samples of J0rgensen et al. (1995). 

For man y of the galaxies, th e raw kinematic profiles reach 
2-3 (see Kelson et al. 1999a). However, due to the effects of 



seeing and the large (metric) slit width, sophisticated modeling 
is required to derive the true kinematic profiles to such large 
radii. 

Eh- As appear to be largely supported by rotation. This is con- 
sistent with the presence of disks, as shown in HST imaging of 
i ntermediate redshift E-hA galaxies in this, and other cluster s 
(|Franx 1993|, |Wirth, Koo, & Kron 1994 |Kelson et al. 1999a|). 



These galaxies do not appear to be massive, having velocity dis- 
persions of (T i 100 kms"^ A more detailed discussion of these 
galaxies is given in the context of the fundamental plane and 
the M/L ratios of cluster galaxies in Kelson et al. (1999b). 

This sample, with more than fifty galaxies in a single clus- 
ter, shows that the new, large-aperture telescopes can provide 
extremely accurate, high quality spectroscopic data on distant 
galaxies. Such data provide estimates of mass scales, such as 
those presented here, and absorption line stren gths, to be dis- 
cussed in a future paper (Kelson et al. 1999c). In the future. 



larger surveys will provide a more coherent picture in which 
several hundreds of such galaxies can be analyzed simultane- 
ously, leading to precision measurements of the evolving mass 
function, constraining the detailed evolution of M/L ratios and 
merger rates in both the field and in clusters. 
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Fig. 1. — The color-magnitude diagram for the cluster in S and Vsoe" ^8 14 (van Dokkum e( fl/. 19980- All 194 confirmed cluster members in the HST mosaic are 
shown. The filled circles represent those galaxies in the fundamental plane sample. 
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Balmer line residuals. Signal-to-Noise estimates are listed in Tablet 



Kelson et al. 



600 
^ 400 

3 

3 200 


1000 



S 500 



212 








1 




6000 6500 



400 



g 200 



1 1 
- 215 








1 



6000 



6500 



- 233 




1,1- 







1000 



500 



234 



6000 



6500 



6000 



6500 



400 

o 200 
J 



1500 
Ml 000 



° 500 



- 236 


1 ' ' 


- 1 1 1 1 


1 


1 , , - 



6000 



6500 



- 



: 256 




1 ' ' J 






1 , , - 



600 

m 

S 400 
" 200 





1500 

i°iooo 
a 

3 

3 500 





6000 



6500 



269 




6000 



6500 



6000 



6500 




- 292 










1 



6000 



6500 



A (A) 



6000 6500 
A (A) 



1000 



- ' 
: 298 












% 




6000 6500 



500 



1 1 
- 300 










1,1. 



6000 



6500 



600 

m 

fi 400 
u 200 




600 

400 
P 

" 200 




303 










1,1- 



500 



309 


1 ' ' 








1,1- 



6000 



6500 



6000 



6500 



328 





400 



p 200 

o 

u 



6000 



6500 



400 


- 343 






5 

g 200 

o 






1 








6000 


6500 



335 


1 ' ' 








6000 6500 



r 1000 



500 




6000 



6500 



d 500 
3 
o 
u 





■ 356 










1 



800 
600 



n 400 
o 

" 200 



359 




1 1 1 _ 






1 



6000 6500 
A (A) 



6000 6500 
A (A) 



Kelson et al. 



17 




300 

200 

g 100 
u 





600 



iS 400 
cS 200 



1000 



§ 500 
o 



6000 6500 


■ 368 












6000 6500 


_ 371 


• 


III- 


6000 6500 


: 375 












6000 6500 



■ 366 


II' 






1 


1,1- 



400 
300 

m 

S 200 
" 100 




400 



200 



- 391 


1 1 ' 








1 , , - 



600 
-2 400 


3 

c3 200 


400 
S 200 



6000 6500 


: 369 
»^ 












6000 6500 


372 






1 1 




Ill- 


6000 6500 


- 381 


1 






6000 6500 



- 397 






1 


1 



6000 6500 
A (A) 



6000 6500 
A (A) 



Fig. 2. (continued) — 
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Fitting method of Franx et al. (1989). The comparison shown was made using the adopted template HD7S24, and uniform pixel weighting. As in Fig. nl the 
median offset is < 1%. We conclude that velocity dispersion derived using the techniques in § p.2| can be compare directly to measurements made using the Fourier 
Fitting method on nearby galaxies. 
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Fig. 7. — (a ^ A comparison of velocity dispersions derived using different template stars. These measurements were made using the direct fitting algorithm 
outlined in 95.3. The dispersions derived from each template are compared here to the average of the dispersions derived from the five G9-K3 templates, (b) The 
reduced ot the fit for each template as a function of velocity dispersion. For galaxy #360, does not vary greatly with template because of poor S/N. 
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Fig. 8. — A comparison of rr measured from the left and right halves of each spectrum. The differi 
of the galaxy spectra are shown as a function of a for a given template. The direct fitting method of §| 
necessary. 



in dispersions measured from the left and right halves 
was used, with pixel-weighting and masking wherever 
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Fig. 9. — Spatial profiles of V, and a. These have not been corrected for seeing effects, or for inclination. 
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Fig. 9. (continued) — 
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Fir, 1 — The Inirarithm ir gradient of the second moment of the kinematics, fi2 = V a-^ + V-^, of nearby galaxies using long-slit observations from the literature 
([Simien & Prugniel 1997alb,c). A first-order least-squares fit to these data is shown as the solid line. A second order least-squares fit to the rms scattpL.is used to 



show the ±1(T and ±2a variation with radius along galaxy major-axis (the dashed lines). The profiles of the CL1358-I-62 galaxies, shown in Figure [ll|, are very 
similar to those shown here of nearby galaxies. 
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Fig. 1 1 . — The logarithmic gradients of the second moment of the kinematics, /t2 = V(t^ as fnnr tinrs of radius in arp pr (a,c), and radius with respect to 
the half-light radii of the galaxy (b,d; from bulge-plus-disk profile fitting to HST surface photometry; see [Kelson et al. 1999a). These have not been coirected for 



seeing effects. Artificial scatter has been added to the radii in (a,c) to aid in displaying the data. The E/SO galaxies are shown ni the top panels, aniUJie spirals in the 
bottom. Within a half-light radius, the kinematic profiles appeal' to be independent of moiphology. The least-squares fit to the local data of Fig. are shown here 
in (b,d) and indicate that the kinematic gradients of the CL1358-I-62 galaxies are similar to nearby galaxies. 
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Fig. 13. — (a) Q, a quality of fit defined as the ratio of galaxy flux to the standard deviation of the residuals of the velocity dispersion fit, as a function of the 
expected S/N per pixel, (b) the ratio of Q/(S/N) as a function of S/N. Note that the higher the S/N, the more the residuals from the velocity dispersion fit become 
dominated by sources of error other than the expected noise. In (c), we plot the ratio of velocity dispersion to its formal error, a/Sa, as a function of Q, and in (d) 
and (e), we plot the ratio of cr/Sa to Q as functions of Q and 7, respectively. The formal uncertainties in the velocity dispersions are reasonably well given by the 
quality of the velocity dispersion fit, though 7, the line-strength, is clearly a third parameter, indicating that template mismatch is playing a role. The five data points 
in the lower left corner of (e) are the E+A and emission line galaxies, which are poorly fit by the G9III template star, even when the Balmer lines are excluded from 
the fit. In these cases, the formal errors of the velocity dispersions are higher than one would expect from the quality of the fit, presumably because of mismatch 
with the late-type stellar template. 
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Table 1 
Stellar Templates 



ID 


Spectral 


Date of 




Type 


Observation 


The Sun (twilight) 


G2V 


5/96 


HD5268 


G5IV 


8/96 


HD 102494 


G91V 


5/96 


HD6833 


G9111 


8/96 


HD72324'' 


G9III 


5/96 


HD56224 


Kim 


5/96 


HD4388 


K3in 


8/96 



Note. — Best-fit template 



Table 2 

Features used to Measure the Instrumental Resolution for the Galaxy 

Spectra 



A 


Species 


A 


Species 


5224.137 


9-2Pl(2.5) 


6533.044 


6-1 Pl(2.5) 


5238.747 


9-2Pl(3.5) 


6553.617 


6-1 Pl(3.5) 


5577.338 


[OI] 


6577.183 


6-1 PI (4.5) 


5589.128 


7-1 Pl(2.5) 


6596.643 


6-1 P2(4.5) 


5605.143 


7-1 Pl(3.5) 


6603.990 


6-1 Pl(5.5) 


5915.301 


8-2Pl(2.5) 


6889.288 


7-2 P2(1.5) 


5924.723 


8-2 P2(2.5) 


6900.833 


7-2Pl(2.5) 


5932.862 


8-2Pl(3.5) 


6912.623 


7-2 P2(2.5) 


5953.420 


8-2Pl(4.5) 


6923.220 


7-2Pl(3.5) 


5970.278 


8-2 P2(4.5) 


6939.521 


7-2 P2(3.5) 


6192.937 


5-0 P2(1.5) 


6948.936 


7-2Pl(4.5) 


6202.720 


5-0Pl(2.5) 


6969.930 


7-2 P2(4.5) 


6221.771 


5-0Pl(3.5) 


6978.258 


7-2Pl(5.5) 


6287.434 


9-3Pl(2.5) 


7003.858 


7-2 P2(5.5) 


6300.304 


[0 1] 


7011.202 


7-2 PI (6.5) 


6363.780 


[OI] 







Note. — 



Wavelengths and identifications taken from Osterbrock et al. (1997). 
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Table 3 
Final Velocity Dispersions 



Adopted 

ID Moiphology S/N per A o"df.f* (km/s) o"df.k* (km/s) 



95 


SO 


48 


219.1 


±6.4 


220.5 


±6.4 


110 


so/a 


71 


171.8 


3.7 


174.8 


3.7 


129 


so/a 


65 


167.0 


3.6 


167.5 


3.8 


135 


SO 


40 


161.2 


5.4 


159.7 


5.1 


142 


so/a 


43 


148.4 


4.5 


147.9 


4.4 


164 


so/a 


51 


181.0 


4.5 


180.6 


4.3 


182 


SO 


28 


132.5 


5.6 


130.1 


5.4 


209" 


Sa 


35 


66.3 


7.7 


73.5 


7.5 


209" 


Sa 


63 


109.0 


5.8 


107.9 


5.2 


211 


SO 


44 


189.1 


6.2 


190.5 


6.4 


212 


e 


53 


182.3 


4.9 


181.0 


5.7 


215 


SO 


33 


162.7 


6.1 


166.2 


5.9 


233 


E/SO 


74 


238.0 


4.3 


234.9 


4.5 


234 


Sb 


42 


113.1 


8.5 


108.9 


8.2 


236 


SO 


43 


197.9 


5.7 


196.7 


5.6 


242 


E 


58 


207.6 


5.2 


212.4 


5.3 


256 


E 


93 


257.8 


4.1 


261.8 


4.3 


269 


E/SO 


92 


317.2 


5.4 


319.5 


5.8 


283 


Un 


34 


147.2 


5.2 


145.7 


5.2 


292 


SO 


56 


111.3 


3.5 


112.3 


3.2 


298 


SO 


85 


296.1 


5.5 


296.6 


5.2 


300 


SO 


60 


254.4 


5.7 


248.2 


5.5 


303 


E 


52 


164.5 


4.0 


163.5 


4.3 


309 


E/SO 


61 


229.9 


4.7 


228.7 


4.7 


328 


Sa 


59 


75.8 


4.7 


81.2 


4.2 


335 


so/a 


37 


128.0 


4.6 


129.1 


4.3 


343 


SO 


36 


63.1 


5.6 


65.1 


5.0 


353 


E/SO 


77 


223.6 


4.3 


228.5 


4.4 


356 


Sa 


64 


191.5 


3.8 


190.2 


3.7 


359 


SO 


56 


198.7 


4.4 


200.3 


4.5 


360 


E 


9 


173.7 


17.7 


177.5 


19.7 


366 


so/a 


43 


226.3 


6.6 


228.9 


6.4 


368 


Sab 


21 


148.2 


7.4 


145.6 


7.5 


369 


so/a 


36 


129.7 


6.3 


128.2 


6.1 


371 


Sa 


59 


175.9 


3.9 


174.8 


4.0 


372 


Sa 


45 


144.7 


5.0 


142.3 


5.0 


375 


E 


78 


307.0 


6.9 


311.3 


7.0 


381 


E/SO 


48 


210.6 


5.6 


212.4 


5.5 


391 


E/SO 


62 


252.9 


7.1 


260.1 


6.3 


397 


so/a 


38 


135.1 


5.2 


135.9 


5.1 


408 


SO 


47 


233.3 


6.2 


237.1 


6.2 


409 


E 


21 


106.9 


6.9 


107.3 


6.7 


410 


SO 


32 


150.9 


5.0 


148.6 


5.2 


412 


E 


41 


170.7 


6.3 


169.4 


6.3 


433 


so/a 


37 


104.4 


3.7 


106.6 


4.0 


440 


so/a 


25 


61.0 


4.8 


63.7 


5.0 


454 


so/a 


63 


149.1 


3.8 


149.6 


3.7 


463 


SO 


57 


284.6 


6.4 


284.4 


6.4 


465 


Sa 


37 


156.3 


4.7 


156.7 


5.0 


481 


SO 


25 


106.8 


6.0 


107.2 


5.5 


493" 


E/SO 


18 


118.7 


9.0 


116.5 


8.3 


493" 


E/SO 


17 


125.1 


10.1 


127.5 


8.8 


523 


so/a 


42 


171.8 


5.9 


174.6 


5.9 


531 


E 


90 


285.3 


4.8 


281.3 


5.0 


534 


E 


30 


120.3 


6.2 


120.9 


6.1 


536 


E 


72 


252.2 


5.0 


254.1 


5.0 


549 


Sab 


30 


79.9 


4.7 


79.9 


4.6 



Note. — Notes: o-qf.f refers to measurements made using an analog of the Fourier 
Fitting method except that the x~ summation is performed in the reaLdomain. (tdf.k 
refers to measurements made using the direct fitting method of § p.2| . (a) Galaxies 
#209 and #493 were observed in two different slit-masks and positi on a ngles The 
first listed value for #209 was derived from a mis-aligned slitlet (see §Ra). (b) These 
values have been corrected for aperture size to a nominal aperture with diameter of 
3".4 at the distance of Coma. 



